Investigation of the role of shell structure in quasi-fission mass distributions 
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Systematic measurements of mass-ratio distributions for fission following collisions of 48 Ti pro- 
jectiles with even-even target nuclei from 144 Sm to 208 Pb have been made at sub-barrier energies. 
They show the presence of quasifission, and depend strongly on target nucleus deformation and the 
fissility of the composite nucleus. A new framework to analyse systematic mass-ratio measurements 
allows direct comparison with the trends expected from shell structure, independent of assumptions 
or fits. This indicates that quasi-fission mass distributions show trends consistent with low energy 
mass-asymmetric fission of the same actinide elements. 
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The elements found naturally on Earth were formed 
several billion years ago, many in violent cosmic events 
such as super novae. The heaviest of these, Thorium 
(atomic number Z=90) and Uranium (Z=92), have un- 
dergone significant decay since their formation. Their 
continuing a-decay provides much of the energy of vol- 
canos, and the motion of the continents over the surface 
of the Earth From energetics, the ground-states of all 
atomic nuclei heavier than the Iron/Nickel elements are 
in principle unstable to nuclear decay. This can be un- 
derstood from the energy change occurring if they were 
formed by fusion of two lighter nuclei - the process costs 
energy, thus the equivalent decay would release energy. 
However, hardly any naturally occurring isotopes have 
measurable decay rates. This is because the probability 
of quantum tunneling through the potential barrier be- 
tween their compact ground-states and lower energy con- 
figurations is negligible, particularly for decay into two 
similar mass nuclei, a process called fission. However, fis- 
sion in principle provides a definite limit to the existence 
of the chemical elements. For /3-stable nuclei, the liquid 
drop model predicts that nuclei with Z>115 will have no 
fission barrier. However, even heavier elements have been 
formed in the laboratory [2[, with lifetimes >ms. 

This is a result of nuclei with close to integer axis ratios 
exhibiting "magic" numbers of neutrons and/or protons 
(closed shells) which have a lower energy than neighbor- 
ing configurations. This increases the effective height of 
the fission barrier, thus providing stability beyond the 
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liquid drop limit. The shell-stabilization of super-heavy 
elements was already calculated in 1966 |3[ to occur at 
neutron number N=184 and Z=114 (more recent calcu- 
lations [4[ predict stabilization also for Z=120, or even 
Z=126). That prediction initiated major experimental 
efforts to reach this "holy grail", to allow testing of nu- 
clear structure models and relativistic effects in chem- 
istry [5]. Elements up to Z=118 have now been cre- 
ated [2], but forming even heavier isotopes and elements 
will be more difficult, despite some being predicted to 
have even longer lifetimes. 

The difficulty is that super-heavy elements are formed 
in collisions of two massive nuclei, whose large Coulomb 
repulsion inhibits fusion. After contact, the dinuclear 
system is pictured as diffusing in shape [6] over the po- 
tential energy surface (PES). The system may reach com- 
pact shapes inside the barrier - defined as fusion. It is 
generally more likely that following contact it quickly 
re-separates into fragments with masses between those 
of the projectile and target nuclei - a process known as 
quasi-fission @, H| - which inhibits fusion. 

The effect of shell structures is central both to the for- 
mation of super-heavy elements, as well as to their ex- 
istence and properties. Theoretical ideas [9[ and exper- 
iments [10[ suggest that collisions of doubly closed-shell 
nuclei give less quasi-fission, thus forming heavy elements 
with the highest probability. Indeed, the heaviest ele- 
ments have been formed [2[ in reactions with the doubly 
closed-shell nucleus 48 Ca. However, the heaviest practi- 
cal target nucleus is Cf (Z=98), thus for a 48 Ca projectile, 
Z is limited to 118. Shell valleys have also been sug- 
gested as a mechanism to form neutron-rich super-heavy 
isotopes through inverse quasi-fission [12| in collisions of 
two deformed actinide nuclei [13]. Shell structures en- 
countered durin g th e diffusion are also expected to affect 
the dynamics [6|, [ll| . The resulting valleys in the PES are 
thought to intercept flux [12] that might otherwise lead 
to fusion, resulting instead in quasi-fission. The role of 
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FIG. 1: (Color online) Representative mass- angle distribu- 
tions (above) and projected mass-ratio distributions (below) 
for the target nuclei indicated. Each projected distribution is 
normalized to 1050 at the maximum, and is offset from the 
lower curve by 500, or 800 for 208 Pb. 



shell structures in fusion dynamics is the most complex 
and open question, and is a key element in evaluating 
prospects [12] of forming more superheavy isotopes, pos- 
sibly using neutron-rich radioactive beams. 

Experimental determination of the effects of shell 
structure, through measuring the characteristics of quasi- 
fission, is made challenging by the many factors affecting 
reaction dynamics. Even assigning the important neu- 
tron and proton numbers that may control quasi-fission 
mass distributions has up to now involved assumptions. 
In this work we propose a new systematic approach. This 
combines extensive measurements of quasi-fission in re- 
actions forming actinide elements, with empirical knowl- 
edge of the neutron and proton numbers controlling low 
energy fission mass-splits of the same elements. It makes 
use of a key feature expected of shell-induced structure 
in the PES: it will be present for many fissioning nuclei, 
at mass-splits where the same fragment is produced. 

To investigate experimentally the role of shells in quasi- 
fission, reactions must be chosen that give a substantial 
quasi-fission yield. In general, quasi-fission is expected 
to be more likely the larger the charge product of the 
two colliding nuclei, but details of the reaction are also 
important. A heavy reaction partner with a large static 
deformation results in a range elongations at contact, and 
thus of entry points into the PES, depending on the ori- 
entation. If the long axis is aligned with the projectile 
nucleus, the dinucleus is elongated at contact, whereas 
it is more compact if anti-aligned. The aligned orien- 
tation can be selected by choosing a sub-barrier beam 




FIG. 2: (Color online) The width (<jmr) of the Gaussian best 
fit to the measured mass-ratio distributions, as a function of 
the compound nucleus charge Z and the quadrupole coupling 
parameter f3 2 (see text). The width for Z=94 ( 178 Hf target) 
is infinite, as the measured distribution is basically flat. 



energy [181-120]. where measurements indeed show the 
broadest mass distributions, which is strong evidence 
for quasi-fission J2CH26|. Considering also the expected 
damping of shell effects with increasing E^, sub-barrier 
energies should be optimal for an experimental investi- 
gation. 

To carry out the measurements, pulsed beams 
(~ 1.5 ns FWHM) of 48 Ti in the energy range 198 - 
235 MeV were provided by the ANU 14UD electrostatic 
accelerator. Measurements were made at sub-barrier en- 
ergies, 0.98±0.01 of the barriers predicted by Swiate- 
cki's model [27j . This both increases the probability of 
quasi-fission and minimizes E^, emphasizing the influ- 
ence of shell valleys in the PES. The beams bombarded 
enriched targets of neutron-rich isotopes of all even-Z tar- 
gets from Sm to Pb except Gd; 144 Sm was also included. 
The targets were ~50 /ig/cm 2 in thickness, evaporated 
onto ~15 jWg/cm 2 C backings (facing downstream), their 
normal angled at 60° to the beam axis. Binary events 
were measured in two 28 cm x 36 cm position sensitive 
multi-wire proportional counters (2l| . Located on oppo- 
site sides of the beam, their scattering angle coverage was 
5° - 80° and 50° - 125°, allowing full efficiency detection 
of all mass-splits between the projectile and target, for 
centre-of-mass angles (# c . m .) between 40° and 140°. 

Mass-angle distributions (MAD) were extracted as de- 
scribed in Refs.[24]-[26j by determination of the velocity 
vector of each coincident particle. The fission mass-ratio 
Mr - the mass of one fragment divided by the total mass 
- was determined event-by-event from the ratio of the 
two fragment velocities in the center-of-mass frame (2lj . 
Representative MAD are shown in the upper panels of 
Fig {TJ The intense bands at Mr ~ 0.2 and 0.8 corre- 
spond to elastic and other peripheral collisions, with the 
fission and quasi-fission events lying between. 

Total mass-ratio distributions were obtained by pro- 
jecting the data for 40° < # c . m . < 140°. Six representa- 
tive distributions of the ten measured are shown in the 
lower panel. At Mr values more asymmetric than those 
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FIG. 3: (Color online) 3-D isometric and 2-D contour representation of the systematics of the experimental mass-ratio (M^) 
distributions, as a function of the atomic number Z of the composite nucleus, for the reactions indicated. The target nucleus @2 
values are plotted for each reaction. The overlaid joined circles and squares on the 2-D map show the positions of the empirical 
neutron and proton numbers correlated with the low energy asymmetric fission of the same actinide elements. 



shown, yields rise rapidly to the elastic peak. Although 
some distributions show structure, to present the gross 
trends they have been fitted in the range 0.35 < Mr < 
0.65 with a Gaussian function. The resulting cjmr values 
are shown in an isometric histogram in FigJSJ as a func- 
tion of both the compound nucleus Z and the absolute 
value of fa coupling the ground-state with the first 2 + 
state of the target nucleus. For ^2 > 0.1, it is reasonable 
to associate fa with the static deformation. For compari- 
son, the a mr values for 16 0-induced fusion- fission (where 
the effect of shells is generally small) range from ~0.06 
to 0.08 for Z=100 |l[. 

The shape of the wide mass distributions seen for most 
48 Ti reactions must be controlled by two variables: (i) 
the probability of quasi- fission and (ii) the mass (Mr) 
distribution of the quasi-fission. The narrowest distri- 
butions for the 48 Ti reactions (see Fig. 3) result from re- 
actions with the light 144 ' 154 Sm nuclei (forming Z=84), 
and with the heavy doubly magic 208 Pb nucleus (forming 
Z=104). Reactions with Ca projectiles on these targets 
have previously been associated with small probabilities 
of quasi-fission (28[. However, quasi-fission is still present 
for 48 Ti+ 208 Pb, as seen in the MAD (FigHJ by the cor- 
relation of fission mass with angle [29|. The <jmr val- 
ues (Fig|2|) for reactions with the lighter well-deformed 
(fa^0.33) prolate target nuclei show a rapid increase in 
o~mr with Z, and thus fissility. However as Z increases 
beyond 94, an irregular decrease in <jmr occurs, as fa re- 
duces. How these trends correlate with the quasi-fission 
probability clearly depends on the quasi-fission Mr dis- 
tributions, which may be affected by shell structure. 

To investigate the role of shell structure, a new rep- 
resentation of systematic experimental Mr distributions 
is introduced, which highlights persistent shell features, 
which can be difficult to isolate in limited statistics sub- 
barrier measurements. Fig |3] shows the experimental dis- 



tributions of fission yield (color scale) as a function of Mr 
(X-axis) and atomic number Z of the composite system 
(Y-axis), for all the measurements made in this work. 
Contour graphs are presented both as a 3-D isometric 
view (left) and as a 2-D map (right). As in FigJTJ the 
yield is normalized to 1050 at the highest point in the dis- 
tribution, to allow easy visualization of changes of shape 
with Z. From 154 Sm to 174 Yb (fa -0.33) the distribu- 
tions start to show mass-asymmetric shoulders (see also 
FigH]), whose yield increases with system fissility. These 
are qualitatively consistent with results for 48 Ca+ 168 Er, 
attributed [33] to quasi-fission. The mass distribution 
is broadest for 178 Hf (fa = 0.28) even showing indica- 
tions of a dip at symmetry. However, for still heavier 
targets, having decreasing fa but increasing fissility, the 
shoulders appear to become narrower. Quasi-fission may 
be decreasing in probability, and/or shell effects may be 
modifying the quasi-fission mass distributions. 

To investigate whether this behavior could be consis- 
tent with the effect of shell structures, we treat within 
the same framework the systematics for spontaneous and 
low-energy fission of isotopes of the same actinide ele- 
ments - which exhibit shell-driven mass-asymmetric fis- 
sion, with different mass-split modes [14|,LL5j. Fig0]shows 
Mr as a function of the atomic number Z of the fission- 
ing nucleus. The large symbols represent the centroids 
of the empirically determined fission modes for the in- 
dicated isotope of each element, taken from Appendix 
A of Ref. [16], and Ref. 0. The "Standard II" fission 
mode (shown by the green squares labeled S II) gener- 
ally has the highest yield. The yellow circles (S I) repre- 
sent the "Standard I" fission mode, except at Fermium 
(Z=100), where it represents symmetric fission found for 
mass number > 258, associated with two fragments close 
to the doubly-magic 132 Sn. Assuming the N/Z ratios 
of the fragments are the same as that of the fissioning 
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FIG. 4: (Color online) Representation of the systematic 
trends of spontaneous and low energy fission mass-ratios for 
the indicated isotopes of actinide elements, as a function of 
their atomic number Z. Large circles and squares represent 
experimental fission modes, whilst smaller symbols represent 
trends for fixed neutron and proton numbers in the fission 
fragments (see text). 



nucleus (N°/Z°), for a given particle number Z Shel1 or 
^Sheii ^ n ^ e f ra g me nt, the associated mass ratio M^ 1611 
is given by the ratio Z Shelt /Z° or N 5/ie "/N°. The joined 
small circles indicate the M^ 1311 values associated with 
the spherical closed shells Z=50 and N=82. The small 
squares show the expected trends associated with Z=55 
and N=86, the proton and neutron numbers empirically 
found to be closely associated [16| with the generally pre- 
dominant "Standard II" mode. The systematic behavior 
of the mass-splits is consistent with the trends expected if 
fixed proton and/or neutron numbers in the nascent frag- 
ments [16[ are responsible. It appears that the spherical 



shells do not play the most significant role in low energy 
mass asymmetric fission. 

To compare with the quasi-fission data, overlaid on 
the 2-D contour map in Figj3]are plotted the Mf> hel1 val- 
ues for the same shell numbers. For the measured dis- 
tributions with Z>94, it is striking that areas of high 
yield away from symmetry show a systematic correla- 
tion with the trends of the shell structures. This sug- 
gests that shell structure in the PES does contribute to 
the observed mass distributions in these quasifission reac- 
tions, although not necessarily the spherical shells. Shell 
structure in both N and Z at the same mass-ratio should 
affect mass distributions most strongly. In these mea- 
surements the empirical Z=55 and N=86 lines are not 
as close together as in low energy fission, because the 
isotopes formed following capture of 48 Ti projectiles are 
less neutron-rich. Their occurrence at different mass- 
splits means that we should not necessarily expect to see 
in these quasi-fission measurements asymmetric fission 
peaks identical to those in low energy fission of more 
neutron-rich isotopes. 

This new approach allows investigation and tracking of 
neutron and proton shell effects in both low energy fis- 
sion and quasi-fission. How the effects of shells reinforce 
each other could be investigated through a systematic 
study of sub-barrier quasi-fission, by varying the com- 
posite system neutron number for fixed Z. If shells in 
both Z and N occur at the same M#, this should result 
in more distinct peaks in the quasi-fission mass distri- 
butions. Multi-dimensional PES calculations will help 
to predict optimal reactions for measurements. Such a 
study would be a good candidate for an early application 
of radioactive beams to better understanding superheavy 
element formation. 
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